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Fumigation  is  generally  accomplished  by  introducing  a  volatile 
liquid  or  a  gas  into  an  enclosure  containing  the  materials  to  be 
fumigated.  Its  effectiveness  depends  on  the  exposure  time  and  the 
concentration  at  a  given  temperature.   Since  fumigation  chambers  are 
not  usually  designed  to  hold  pressures  above  atmospheric,  the  vaporizing 
liquid  or  the  compressed  gas  that  is  introduced  displaces  an  equivalent 
volume  of  air-fumigemt  mixture.  Losses  of  fumigant  due  to  other  factors 
such  as  leakage  or  wall  absorption  also  occur,  and  if  these  are  known, 
the  concentration  of  the  fumigant  at  any  time  can  be  calculated. 

This  discussion  applies  only  to  chambers  operating  under  atmospheric 
conditions,  which  contain  only  air  initially  and  remain  at  atmospheric 
pressure  throughout  the  fumigation  period.   Vacuum  fumigation  will  not 
be  considered. 

The  concentration  of  the  fiomigant  in  the  chamber  at  any  time  is  a 
function  of  two  factors — the  quantity  introduced  per  unit  volume  of 
chajnber  space,  and  the  rate  of  leakage  through  the  walls  and  crevices 
of  the  chamber,  both  as  displacement  loss  during  input  and  by  diffusion 
to  the  outside  atmosphere  during  the  fumigation  period.  As  a  basis  for 
calculation  it  is  assumed  that  the  following  conditions  are  maintained: 
(1)  the  rate  of  input  of  the  fumigant  is  constant,  (2)  the  input  time  is 
much  shorter  than  the  entire  fumigation  period,  (3)  the  fumigation 
mixture  in  the  chamber  is  in  the  gas  phase  and  homogeneous  at  all  times, 
(4)  the  components  of  the  mixture  are  always  in  the  same  proportion 
(except  as  will  be  mentioned  later),  (5)  the  walls  of  the  chamber  are 
rigid  and  allow  slight  interchange  between  the  gases  in  the  chamber  and 
the  outside  atmosphere,  and  (6)  the  leakage  rate  of  gas  (fumigant  plus 
air)  to  the  outside  atmosphere  is  constant. 
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As  fumigant  is  discharged  into  the  chamber,  it  vaporizes  and  mixes 
with  the  air  in  the  chamber,  displacing  an  equal  volume  of  air  to  the 
outside  atmosphere.  The  gas  displaced  initially  consists  only  of  air, 
but  it  gradually  becomes  richer  in  fumigant,  until  the  concentration  of 
incoming  fumigant  is  in  equilibrium  with  that  of  the  gas  being  displaced. 
The  concentration  upon  completion  of  input  will  thus  be  a  function  of  the 
volume  and  rate  of  fumigant  input  and  the  quantity  lost  by  displacement. 

After  the  input  of  fumigant  is  completed,  there  is  no  longer  any 
displacement  loss.  Since  the  fumigation  period  is  much  longer  than  the 
input  time,  leakage  loss  is  usually  the  more  important  factor  at  this 
period.  Rate  of  leakage  of  fumigant  is  highest  at  the  time  of  its 
maximum  concentration,  and  then  gradually  diminishes.  Since  the  chamber 
gases  are  in  continuous  motion  as  a  result  of  fan  circulation,  the 
leakage  rates  of  different  gases  are  more  or  less  the  same.  However, 
when  differences  in  leakage  rates  of  the  components  of  the  fumigant  are 
observed,  the  gas  whose  leakage  rate  appears  to  be  abnormally  high  may 
be  adsorbed  on  the  interior  surfaces  or  removed  in  some  other  manner, 
such  as  chemical  reaction  with  materials  in  the  chamber,  decomposition, 
or  condensation. 


Equations:^ 

Based  on  these  considerations  and  by  means  of  the  calculus,  equations 
have  been  developed  to  express  the  concentration  relationships.   The 
symbol  notations  and  mathematical  derivations  of  these  equations  will  be 
found  in  the  appendix  of  this  report.   The  equations  are  similar  to 
those  of  monomolecular  or  monoatomic  reactions,  such  as  radioactive  decay, 
in  which  the  reaction  rate  (leakage)  is  proportional  to  the  concentration 
of  reactant  (fumigant)  at  any  moment  during  the  reaction  (fumigation). 

During.  Fumigant  Input 

5i  =  (100/Z)(1  -  e-Z^)  +  Soe-Ze  (l) 

e  =  (2.303/z)  log  /Tioo  -  zSo)/(ioo  -  zsiX^   (2) 

Usually  Sq  =  0,  and  the  equations  become 

Si  =  (100/Z)(1  -  e-2e)  (IT) 

e  =  (2.303/z)  log  ^00/(100  -  zsi^        (2O 

Fumigation  Period 

log  S  =  log  Si  -  0.4343  kt»/V  (3) 


ij   The  ten  text  equations  (1,  2,  1»,  2»,  3,  4,  5,  6,  7,  and  8)  are 
derived  in  the  appendix,  where  they  are  designated  A-I4,  A-I3,  A-14S 
A-13S  A-I7,  A-34,  A-31,  A-22,  A-24,  and  A-26,  respectively. 
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Interpretation  and  Utilization  of  Equations 


Quantity  of  Fumi^ant  Dispensed  and  Initial  Concentration.   From 
equations  1  and  2  the  quantity  of  fumigant  dispensed,  Q,  and  the  initial 
concentration.  Si,  can  be  calculated  from  each  other.  The  volume  of 
fumigant  gas  input  can  be  estimated  from  the  weight  of  fumigant  dis- 
charged and  the  ideal  gas  law.   Input  rate,  I,  can  be  determined  by 
timing  the  input.   In  experiments  at  our  laboratory  the  fumigant  was 
dispensed  from  pressure  cylinders,  and  the  cylinder  was  weighed  before 
and  after  discharge  and  the  time  of  input  noted  in  order  to  calculate  I. 
The  leakage  rate,  k,  must  be  determined  separately  for  the  chamber  in 
question;  it  is  best  done  in  advance  by  chemical  analysis.   It  was  found 
convenient  to  discharge  dichlorodifluoromethane  into  the  chamber  and  to 
determine  it  by  means  of  thermal-conductivity  measurements  (l)  at 
frequent  intervals.  Methyl  bromide  (/j.)  and  similar  gases  could  also 
have  been  used,  k  can  be  calculated  from  the  slope  of  a  plot  of  log  S 
against  t*.   If  k  is  small  and  the  input  time  is  short,  k  may  be  neglected 
(Z  =  1)  for  approximating  Q  and  S]_  from  each  other. 

I^feiximum  Obtainable  Concentration  (S^  =  O).  From  equation  1^  we  note 
that  the  maximum  possible  concentration  is  less  than  lOO/Z. 

Input  Rate.   From  equation  1  we  can  estimate  that,  for  a  given  fumi- 
gation dosage,  the  faster  the  input  rate,  I,  the  lov/er  will  be  Z,  and 
consequently  the  greater  will  be  the  concentration  at  the  completion 
of  the  input.  Thus,  to  increase  initial  fumigant  concentration,  the 
input  rate  should  be  as  fast  as  possible. 

Maintenance  of  a  Constant  Fumigant  Concentration.   It  may  be 
necessary  or  desirable  to  maintain  a  fumigant  concentration  at  a  con- 
stant value  after  the  input  is  completed.  This  may  be  accomplished  by 
means  of  a  continuous  secondary  input,  I2,  of  fumigant  into  the  chamber 
to  balance  the  loss  due  to  leakage.  The  rate  of  the  secondary  input  is: 

I2  =  k  Si/(100  -  Si)  iU) 

Half-life  of  Fumigant  Concentration.  The  time  required  for  the 
concentration  of  fumigant  to  fall  to  half  its  original  value  may  be 
calculated  as  follows: 

t»i/2  =  0.693  V/k  (5) 

Plotting  of  Concentration-Time  Curves.   It  is  useful  to  plot  log  S 
against  t,  since  from  equation  3  ^e  can  see  that  the  relationship  will 
be  linear.  The  slope  of  this  line,  mi,  is  equal  to  -  0.4343  k/V,  and 
we  can  therefore  calculate  k  from  mi. 


k  =  2.303  Vmi 
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Additional  Meanings  of  k.  We  have  defined  k  as  the  rate  of  leakage, 
and  have  assumed  that  the  fumigant  is  inert  and  is  not  removed  from  the 
chamber  in  any  other  way.   In  certain  cases,  however,  this  may  not  be  so, 
and  k  includes  all  avenues  of  loss  of  fumigant.  We  may  then  redefine  k 
as  the  sum  of  all  these  loss  rates. 

k=£ki 

where  kj^  =  any  individual  loss  rate,  such  as  that  due  to  leakage,  absorp- 
tion, condensation,  or  chemical  reaction. 

Determination  of  Exposure.  Exposure  (frequently  referred  to  as  Ct) 
is  usually  measured  by  plotting  concentration  against  time  over  the  fumi- 
gation period  and  then  estimating  the  area  under  the  curve.  This  may 
also  be  done  by  cauLculations  based  on  the  equations  that  have  been  derived. 
It  is  convenient  to  split  up  exposiire,  E,  into  two  parts,  Ej  during  input 
and  Ejj  the  fumigation  period.  V/e  shall  consider  only  the  usual  case 
where  the  fumigation  is  started  at  zero  concentration  of  fumigant.  Since 
it  is  customary  to  express  E  as  a  function  of  weight  concentration  of  the 
active  ingredient  per  unit  volume  (milligrams/liter)  and  time  (hoiu?s),  a 
constant.  A,  and  the  factor  I/6O  are  introduced  to  effect  these  conver- 
sions, respectively,  from  percent  concentration  and  time  in  minutes. 

Ei  =  AV  (100  8  -  Si)/60  (I  +  k)  (6) 

Eji  =  AV3]_N/60  k  (y) 

where  N  =  1  -  e'^^VV 

ektVv  =  iog-1  ^o.k3k?>   ktVv) 

and  E  =  Ej  +  Eji 

From  these  equations  we  may  also  calculate  the  necessary  input,  Q,  for 
any  desired  E  at  a  given  time  period,  t».  Usually  Ej  is  smaller  than 
Eji  and  may  be  neglected  for  approximation  purposes. 

.*.  E  =  Eii 

e  =  (2.303/z)  log  ^.67  anv/(i.67  anv  -  kz  Eil17  (8) 

Calculation  of  E  from  Analytical  Data.   If  we  can  neglect  Ej, 
we  merely  have  to  calculate  Eji,  as  follows:  Plot  log  S  against  time, 
and  determine  k  from  the  slope,  ra.  Extrapolate  the  curve  back  to  t»  =  0 
(at  time  of  completion  of  input)  and  estimate  log  S]_  and  hence  S^. 
Knowing  V,  t*,  and  A,  we  can  calculate  Eji  from  equation  7. 

Maximum  Obtainable  E.   From  equation  8  we  can  estimate  the  maximum 
obtainable  E  to  be  1.67  A^rvVkZ. 
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Application  of  Equations  to  Experiment 

In  the  course  of  fumigation  work  in  cooperation  with  the  Army 
Chemical  Corps,  in  tests  of  the  efficiency  of  ethylene  oxide  as  a 
sterilizing  agent,  data  were  collected  which  may  be  used  in  applying 
the  equations  just  derived.  The  data  are  from  two  fumigation  experi- 
ments at  90°  F.  with  60  pounds  of  10  percent  ethylene  oxide  plus  90 
percent  carbon  dioxide  in  the  first,  and  50  pounds  of  12  percent  ethylene 
oxide  plus  88  percent  dichlorodifluoromethane  in  the  other. 

The  carbon  dioxide  and  dichlorodifluoromethane  acted  as  propellants 
to  dispense  the  ethylene  oxide,  and  also  as  flame  suppressants  to  pre- 
vent a  fire  hazard  (2).   From  time  to  time  samples  from  the  gas  mixture 
in  an  approximately  977-cubic-foot  airplane  fuselage  used  as  a  fumigation 
chamber  were  analyzed  chemically  (l).  Results  are  shora  in  table  1. 
Plots  of  the  log  of  concentration  against  fuiriigation  time  are  shown  in 
figure  1. 

Table  1, — Concentrations  (volume  percent)  of  ethylene  oxide  and  propellant 
in  977  cubic-foot  airplane  fuselage.  Temperature  90°  F. 


Time 
(minutes) 


Ethylene  oxide 


Propellant 


Total  furaigant 


Ratio  propellant 
to  ethylene  oxide 


Ethylene  oxide  +  carbon  dioxide  at  input  rate  of  3,874  liters/minute 


oi/ 

02/ 

4.2 

37.6 

41.8 

9.0 

3.6 

35.1 

38.7 

9.8 

7 

3.5 

34.3 

37.8 

9.8 

35 

3.0 

31.3 

34.3 

10.4 

65 

2.4 

27.5 

29.9 

11.4 

152 

1.6 

19.7 

21.3 

12.3 

k 

149 

96 



— 

t]_/2  (min.) 

128 

182 

— 

— 

Ethylene 

oxide  +  dichlorodifluoromethane 

at  input  rate 

of  568 

liters/minute 

t^ 

4.9 

13.2 

18.1 

2.671 

4.6 

14.3 

18.9 

3.1 

15 

4.4 

14.1 

18.5 

3.2 

45 

3.8 

9.9 

13.7 

2.6 

75 

3.0 

10.1 

13.1 

3.4 

165 

2.3 

7.2 

9.5 

3.1 

k 

118 

110 

__ 

— 

^1/2  ("li^*) 

162 

169 

— 

— 

1/  Calculated  from  theoretical  equations. 
2/  By  extrapolation  of  analytical  data. 
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From  the  data  certain  observations  may  be  made: 

Although  equal  weights  of  ethylene  oxide  (6  pounds)  were  used  in 
each  fumigation  in  the  ssime  airplane  fuselage,  the  mixture  containing 
dichlorodifluoromethane  resulted  in  a  higher  concentration  of  ethylene 
oxide.   This  is  in  accordance  with  equation  1,  and  is  due  to  the  greater 
volume  of  total  fumigant  discharged  into  the  airplane  by  the  carbon 
dioxide  mixture.  At  higher  concentrations  this  difference  would  increase 
greatly.  Greater  quantities  of  the  ethylene  oxide-carbon  dioxide  mixture 
than  of  the  ethylene  oxide-dichlorodifluoromethane  mixture  will  be  needed 
to  reach  a  given  concentration  of  ethylene  oxide,  and  the  latter  mixture 
will  also  be  able  to  reach  a  greater  maximum  concentration  in  a  given 
chamber  no  matter  how  much  fumigant  is  dispensed.   The  initial  concen- 
trations calculated  from  the  theoretical  equations  (based  on  an  estimated 
k  value  of  100)  are  in  good  agreement  with  the  results  found. 

The  plots  of  log  concentration  against  time  are  linear,  as  predicted 
by  equation  3« 

The  k  vsilues  for  both  propellants  were  of  the  saime  order  of  magnitude, 
suggesting  that  their  loss  mechanism  was  of  the  same  type.  The  k  values 
for  ethylene  oxide  were  higher,  however,  especially  in  the  first  experi- 
ment, perhaps  because  of  loss  by  absorption  on  the  walls  and  other 
internal  surfaces  of  the  chamber. 
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Notat 

■ion 

Unit 

Q 

liters 

V 

do. 

X 

do. 

L 

do. 

I 

liters/minute 

k 

do. 

t 

minutes 

t» 

do. 

t» 

1/2 

do. 

S 

percent  by  volume 

So 

do. 

Si 

do. 

G 

— 

Z 

— 

E 

mg. -hours/liter 

C 

mg. /liter 

do. 


m 
N 


log  {%   by  vol.)/minute 


.kt»/V 


Appendix 

Meaning 

Volume  of  input  gas 

Voliome  of  fumigation  chamber 

Volume  of  fumigant  gas  in  chamber 

Volume  of  outgoing  gas  due  to  leakage 

Input  rate  of  fumigsmt 

Gas-loss  rate,  k  =  -2.303  Vm 

Time  (0  at  start  of  input) 

Time  of  fumigation  after  input 

Half-life  (S  =  S1/2) 

Concentration  of  fumigant 

S  at  beginning  of  input 

S  at  end  of  input 

e  =  Q/V 

Z  =  (I  +  k)/l 

Exposure  measure,  or  Ct  =  S  x  t  x  A/6O 

Concentration  of  active  ingredient 
(C  =  SA) 

Conversion  factor  from  percent  by 
volume,  A  =  C/S 

Slope  m  =  0.4343  k/V 
N  =  (ektVv  _i)/e-ktVV 
ektVv  =  iog-1  (0.4343  ktVv) 
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Derivation  of  Equations 

During:  Input,   Change  of  volume  of  fumigant  gas  with  time  in  the 
fumigation  chamber: 

^  =  dS  -  x.dQ  _  k.dL  /,  -,  N 

dt   dt   V  dt   V  dt 

where  SS  =  input,  ^»^  =   displacemt  t,  and  iS.2±:  =  leakage 
dt     ^  *  V  dt      ^       *     V  dt 

But  I  =  ^  (A-2) 

dt 

and  K  =  §7^  (A-3) 

dt 

Substituting  for  ^   and  -rr   in  equation  A-1 
dt     Qt 

4r  =  I  -  xl/V  -  xk/V  (A-4) 

dt 

Also  3  =  lOOx/V  (by  definition)  (A-5) 

Differentiating  s  with  respect  to  x 

dx  =  (V/100)  ds  (A-6) 

Substituting  for  x  and  dx  in  equation  A-4 

^  =  100  I/V  -  S  (I  +  k)/V  (A-7) 


and  rearranging  and  setting  up  for  integration 
Integrating  between  limits 


(Si  /tl 

\      ds/^00  I  -   (I  +  k)  S7  =1       (1/V)   dt  (A-8) 


^l-mi^^)7ln^lll:[l:l]l^  (A-9) 

But  Q  =]       I  dt  =  Iti  (A-10) 
•o 

t-i   =  Q/I  (A-11) 
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Substituting  for  t]_  in  equation  A-9  and  converting  to  common  logarithms. 

But  e  =  Z/V  and  Z  =  (I  +  k)/l 
Substituting  these  values  in  equation  A-12 

e  =  (1/Z)  In  ^00  -  Z  So)/(100  -  Z  S]J7  (A-13) 

and  solving  for  S]_ 

Si  =  (100/Z)(1  -  e-2e)  +  s^e-Z©  (A-14) 

If  Sq  =  0,  equations  A-13  and  A-14  become 

9  =  (1/Z)  In  2^00/(100  -  ZSi)y  (A-13') 

and  Si  =  (100/Z)(1  -  e-Z&)  (A-140 

Fnnrigation  Period;  Since  the  input  is  completed, 
1=0,  and  equation  A-7  becomes 

^  =  Sk/V  (A-15) 

dt 

Rearranging  and  setting  up  for  integration  between  limits 

.S  .t» 

\     (1/S)   ds  =  -(k/v)  \     dt  (A-16) 

^Si  -'o 

Integrating  equation  A-16 

log  S  =  log  Si  -  0.4343  ktVV  (A-l?) 

since  k  =  2.303  Vm,   equation  A-l?  becomes 

log  S  =  log  Si  +  mt»  (A-17») 

or  S  =  Sie2.303  mt'  (A-17»0 

Determination  of  E  During  Input.  By  definition  Et  =  area  under  curve 
of  concentration  vs.  time 

/.  Ej  =•(  (C/60)  dt  (A-18) 

Jo 
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But  C  =  AS,  and  substituting  for  C  in  equation  A-18 

Ej  =  (A/60)  1   Sdt  (A-19) 

^o 

From  equation  A-11 

t  =  Q/i  =  ev/i  (A-n) 

Differentiating  equation  A-11  with  respect  to  I, 

dt  =  (V/I)  de  (A-20) 

Substituting  for  S  from  equation  A-li+,  and  t  and  dt  from  equations 
A-11  and  A-20  in  equation  A-19, 

(9 
Ej  =  (AV/60I)\  /Sge-QZ  +  100  (e®^  _  l)/Ze®27  de         (A-21) 

After  rearranging  equation  A-21,  and  integrating  between  limits 
(assuming  Sq  =  O) 

Ej  =  AV  (100  0  -  S;l)/^0  (I  +  k)  (A-22) 

Determination  of  E  During  Fumigation  Period.  Equation  A-l?*'  may 
be  converted  to 

S  =  Si  e-ktVV  (A-17»»») 

Substituting  C  =  SA  for  G  in  equation  A-18, 

-t» 

Eii=(A/60)[     Sie-ktVv  dt  (A-23) 

Integrating  between  limits  and  substituting  N, 

%I  "  ^^^1  ^/^^^  (A-24) 

Also  E  =  Ej  +  Ejj 

Substituting  the  value  for  S]_  from  equation  A-14*  in  equation  A-2A.,  and 
solving  for  6, 

e  =  (2.303/Z)  log  ^.67  ANV/1.67  ANV  -  kZEjj/  (A-26) 

Under  most  conditions  the  input  time  is  short  and  Ej  can  be  neglected. 

E  =  En 
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Maximum  Obtainable  E  for  a  Given  t^.   In  order  for  equation  A-26 
to  have  meaning,  the  quantity  1.6?  AW  -  kZ  Eji  must  be  greater  than  zero. 
The  limit  value  of  Eji  can  be  calculated  by  equating 

1.67  ANV>kZ  Eii  max 

•  Ett  max-:;:  (1.67  MJV/kZ)  (A-27) 

Maximum  Obtainable  Fumigant  Concentration  (S  max).   If  we  assume 
Sq  =  0,  we  may  apply  the  reasoning  just  used  for  equation  A-26  to  equation 
A-13S  and  100 >Z  S]_ 

•  Si  max -<1(100/Z)  (A-28) 

Half-life   (t»  3/2). 
By  definition  S  =  St/2 
Substituting  in  equation  A-17*, 

log  (S1/2)   =  log  Si  +  m  tt^y^  (^-29) 

On  rearrsinging,   this  becomes 


1/2 


m 


t»W2  =  -  0.30103/m  =  0.693  V/k  (A-31) 

Continuous  Secondary  Input  to  Compensate  for  Leakage.  To  accomplish 
this,  equation  A-7  becomes 

ll  =  0  =  (100  I2/V)  -  /si  (I2  +  k)/v7            (A-32) 

Rearranging, 

100  I2/V  =   (I2  +  k)  Si/V  (A-33) 

.*.   I2  =  k  Si/(100  -  Si)  (A-34) 


